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In this work a mathematical apparatus for method of several global Lyapunov functions is applied 
to study the stability properties of nonlinear model of concrete singularly perturbed mechanical 
system "inverted pendulum on a cart" with discontinuous relay- type control. In this analysis, the 
concept of solution in sense of Gelig et al. 8] is used. 

PACS numbers: 



I. INTRODUCTION 



The nonlinear controlled model system " inverted pendulum on a cart" is one of intensively studied last 
time (see, e.g. The results of investigation of this system are applicable at analysis of a number 

of concrete problems of stability, e.g. monorail 0], fault tolerant control inverted pendulum on a 
rotating arm as well as satellites and underwater vehicle with internal rotors etc.pj. 

In the work |{| it was performed a study of stability of nonlinear system "inverted pendulum on a 
controlled cart" for a number of cases. There the conrol was continuous and was taken as sum of the 
controls 



(i) 



and it was applied to the cart moving on the plane surface. 

In recent work of Bloch et al. it was performed the investigation of stabilization of this mechanical 
system using the method of controlled Lagrangian. The control law was obtained in the form 

kbsinf3(a/3 2 + ccos/3) 
" ~ (2) 



a — b 2 /m(l + k)cos[3 

where b = ml, c = —mgl. m is the pendulum mass, I is the pendulum lchgth, g is the acceleration due to 
gravity, m = M + m, M is the cart mass. And equilibrium j3 = /3 = r = r = is achieved if dimensionless 
constant 

am ~b 2 M 

k > — — - — > 3 

b z m 

The control force acts to a cart while there are no direct action to a pendulum. At analysis it is used a 
character of the Lie symmetry group for nonlinear system, which in the case of plane system "inverted 
pendulum on a controlled cart" appears to be that of translation. 

Also in that work it is studied an asymptotic stabilization of the system "inverted pendulum on a 
controlled cart" , and control under investigation is breaked into " conservative" and " dissipative" piece 

u = u cons + u dls , (4) 

which are analysed separately. 

Introduced in the work Lyapunov function for system with control is taken so that its total time 
derivative was non-negative elsewhere (dV/dt < 0) and disappeared at the set M, which is determined 
by equating to zero of dissipative control component u dls = 0. 

As follows from simulation results for concrete parameters of nonlinear system " inverted pendulum on 
a controlled cart" (for given stabilizing action with added dissipation), obtained using the MATLAB 
system, the pendulum begins the motion from almost horizontal position what indicates that attraction 
region is large enough, even at given positive (down) initial velocity. Initial position of the cart is s(0) = 0, 
and initial velocity is s = —3 m/sec. The cart comes to an equilibrium state in the position which is far 
enough from the initial one. The control law was taken with an initial peak to provide the initial large 
acceleration. The Lyapunov function is negative in the beginning and then it strictly increases up to zero 
at the equilibrium. 
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In present work, it is performed a detailed analysis of the behavior of the controlled mechanical system 
"inverted pendulum on a cart" with discontinuous relay- type control on the basis of the method of 
investigation of the global asymptotic stability of nonlinear dynamical systems with using of two Lyapunov 
functions. The distinct feature of such analysis from the case of continuous right-hand sides || is that total 
time derivatives of the Lyapunov functions appears to be also discontinuous. Because of this to study the 
global stability of this system it is used the method of two Lyapunov functions (see 0) which permits to 
perform the stability analysis for discontinuous systems, in particular, for the case of mechanical systems 
with dry friction (see, e.g. 0]). At such analysis the solution is considered in the sense of Gelig et al. 



II. DYNAMICS OF CONTROL OBJECT IN NEW VARIABLES 



It is considered the problem of global stabilization of inverted pendulum on a controlled cart in the pres- 
ence of unknown disturbance acting at cart. This mechanical dynamic system is described by nonlinear 
system of differential equations in the form (see, e.g. Brusin et al. @) 

Mr + L(3cosf3 + Nr - L(3 2 sin(3 = Gu(t) + D{t) , , 

Lcosflr + 1(3 + c(3 + nrcosfl — Lgsin(3 — 0, 

where D{t) is uniformly bounded (D(t) < II where IT is known value), continious function of t, describing 
the external action at cart; rh — M + m, L = ml, c = cq + kI, I = J + ml 2 where M > 0,m > are 
the cart and the pendulum masses, respectively. N > and k > are the coefficients of friction force 
for motion of cart and pendulum, respectively; cq is the coefficient of elastic force moment for rotation 
friction of the pendulum; g is the acceleration due to gravity; G > is the coefficient of amplification 
of motor; J > is the inertion moment for pendulum relative to the mass center; r is the coordinate of 
cart mass center; |/3| < n/2 is the angle between pendulum axis and the vertical, numbered from vertical 
unstable positon of pendulum; u is the value of controlled signal of regulator. 

The purpose of control is to reduce controlled "cart" in asymptotics to given position and attached 
"pendulum" in vertical position (at which its gravity center is situated above fixed point) from any 
initial position (such that \f3\ < ir/2) of "cart" and "reversed pendulum" in the presence of immeasurable 
perturbation at cart, i.e. /3 — > 0, r — > as t — > oo. 

Let introduce new variables (s, s, f2, (l) for state vector of the object of control using relations 

s = r + gQ((3), Sl(/3) = - In | tan(^ - |)|, p G (-|, |), p > 0. (6) 

(3 = 2(tt/4 - arctg(exp(-Q))). (7) 
The variables (s, s) are connected due to system (5) by the following equation 

1 11 D(t) 
(p + a)[^{(3)(s + s) + bs\ =22XiUi + b(s + as)-u (8) 



ip(P) = MI - L 2 cos 2 (3, 

where p = d/dt; a > 0; a > 0; b > are some numbers. A — \G(I — pL)\ > 0; Ui are the continuous 
functions of new variables 



ui = r, «2 = rcos 2 /3, U3 — $ 2 sinf3, 114 = $cosP,u^ — sin2{3, 

u 6 — $/cos{3 = Cl, uj = tg/3, us = as + p/3 2 sin/3 /cos 2 (3, ^ 
ug = (J3sin2f3 — acos 2 /3)s — p/3 2 sinf3, u\q = s + as, 
u n = (p + a){scos 2 (3) = (s + as)cos 2 (3 — ssin2/3. 
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Here functions of old variables are expressed via new ones with taking into account the relations directly 
following from (6) 

r = s- pti(0), f3 = ncos(3, (10) 

and angle (3 is expressed from (7). 
The constants Ai is in the form 



_ pM k + IN- pNL Lk L 

Ai = j ; A 2 = -^-; X ^~q' 

Lc L 2 g pMc pLMg , . 

A 4 = - T ; A 5 = -— ; A6 = -— 5 A 7 = — (11) 

Ml _ L Mia _ ah 1 

As — ~r~ j Ag — — ; Aio — -. — ; An — — — 



Further, if the relations 



A ' ,J " A ' 



r = s — pfl, (3 = ti,cos(3 — ti 2 sin (3 cos (3, 



following from (10), will be taken into account, then equation for Q can be expressed in new variables as 
1 



Ls + KS 



where new designations 



Q + d£l + d 2 tgl3 + d 3 n z sin(3 = (12) 

pL — I 



k P -c Lg I 

dl = ^l> d2 = ^I> d3 = ^J—J (13) 



are introduced. The obtained equations (8) and (12) describe the behavior of initial system in new 
variables. 

Let take the control law in the form [5| 

u(t) = u(t) + Au(t), (14) 

where function [f| 

n 

u(t) = ^\iUi(t) (15) 

i=l 

is taken so that right-hand side of (8) was equal to zero if Au(t) = 0. Au(t) - correcting term to 
compensate perturbation D(t), which is a piecewise function with a jump discontinuity along some 
smooth suface. 

Then (8) for closed system at b = will be wrote in the form (see p|) 

(p + a)[jm(s + as)} = -Au(t) (16) 

Let denote 

^(/3)(s + as)= 7 (t) ; (17) 

then obtain equation for s 

Substitution of (17) in (16) results in 

(p + a )l = -Au(t)-^- (19) 
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7 = -ay - AAu(t) - BD(t), 



(20) 



B = A/G. 

The equations for y — col(Cl, O) (see (12) can be obtained, if express s from (18) 



yip — yip 



as 



or, with taking into account (18) 



yip — yip y 

s = 



ip' 



For numerator of right-hand side of (12) with taking into account (21) one obtains 



Ls + ks = -^(yip — 7V0 + (k — aL)s 



or with taking into account (20) 



Ls + ks = -^((-ay - AAu(t) - BD(t))ip - yip) + s(k - aL) 



Substituting (24) in (12) one obtains 



fl = -diO - d 2 tgp3 - d 3 il 2 sin(3+ 

, 1 , T (-ay-AAu(t)-BD(t))iP~yiP 

(L — h s{k - aL) 



pL-I 



ip 2 



(21) 



(22) 



(23) 



(24) 



(25) 



ip = -L 2 sin2f3. 

So, the dynamics of the control object (5) in new variables s, 7, O, tl will be determined by system of 
four differential equitions of the first order 



' ds 7 



dt 
dfl 
~dt 
dtl 
~dt 

+ 



dt ip(/3) 
dy 



as 



ay - AAu(t) - BD(t) 



n 



= -difl - d 2 tg/3 - d 3 h 2 sinj3+ 

L (-ay - AAu(t) - BD(t))ip - yip 



(26) 



I pL-I 



tl- 2 



+ s(k — aL) 



III. THE CASE OF CONTINUOUSLY ACTING PERTURBATION D(t) 



To stabilize the system it will be used a discontinuous relay-type function 



1, 7 >0 

Au(t) — Tlsigny, signy = { — 1, 7 < (27) 

(-1,1), 7-0 
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So, as a control it will be used a function 

u = u + Au(t), 

where u is taken from (15). 

If denote x = 00/(5,7), y = col(il,Cl), then system (26) can be attributed to the class of systems 
considered in ■ 

To use theorems 1-3 from Q there is necessary the verification of all conditions and proposals used. 
It can be obtained that conditions 1.2, 1.3 of theorem 1 from Q are fulfilled if functions V, W, tp, 77 are 
taken in the form 



f V(s,j) = s 2 + k 1 2 , 



cos?- 2 0(0.) cos?" 1 0(SI) 
^ = £ (|.s| 2 + | 7 | 2 ) 

2afl 2 
V ~ cos?" 2 p(n) 



-r > 



(28) 



where k > 0, e > are the numbers taken by corresponding manner q — 2gL(gL — I) 1 > 1, r 
2(q-l)- 1 Lg( Q L-I)- 1 >0. 



IV. THE CASE OF IMMEASURABLE VELOCITIES: REDUCING TO SINGULARLY 

PERTURBED SYSTEM 

Consider system (5) at D(t) = 0. Let that in the regulator velocities r,B are measured inexactly, and 
their approximate values i\ , z 2 are generated with using the system 0, ll(| . 



~dt 

dz 2 
~dt 



= r- z\, 
= 0-22 



(29) 



where < /i < 1 is the small parameter. The processes approximate processes r,(3. It is essentially 
that in contrast to later they can be calculated without using of differentiation operation 

h(t) = Dl(r{t)), z 2 (t) = £>i (/?(*)). (30) 
In this case instead of globally stabilizing control u = ^ Ajtij(i) consider control in the form 



= XjUjjt) 



(31) 



where Ui is determined from continuous functions Ui by replacement f, to Z\,Z2, respectively. 

U\ = Zi, U2 = zicos 2 f3, 1*3 = z 2 sin0, 114 = z 2 cos(3, 
u 5 = u 5i u 6 = z 2 /cosf3 = ft, M7 = Uf, u s = az 3 + pz\sin(5 j cos 2 (i 

7T 8 

u 9 = (z 2 sin2/3 — acos 0)zs — pz 2 sinf3; Uio — Z3 + ar — apln\tg( — — — |; 

7T i3 

Mil = Z3Cos 2 (3 + sz 2 sin2/3 + ar — apln\tg( — — — )\cos 2 j3. 



Let introduce a vector process z(t) — col(z\(t), Zzit)) 

Zi = zi — f, z 2 



h-0 



(32) 



(33) 



G 



and also denote 



u7 = u. - Ui 



(34) 



Then system (26), (29), (31) with taking into account (33), (34) can be presented in the form of six 
diferential equitions of the first order (cf. with (5)) 

~jT +as = 7/V>(/3) 



d'y 

— + 

dt 

dVL _ 
~dt ~ 
dtl _ 
~d~t ~ 

dz\ 



i=l 



(35) 



dt 



/n(s,7,/z), 



dz 



From above presented it's seen that /n, f r , f/3 are the continuous functions . Moreover, at /x = system 
(35) is reduced to system (26) (at D(t) = and Ati = 0). 

To obtain functions f r and fp in direct form express Zi, i\ z 2 , z 2 via z\ z 2 



z\ = z\ + f, z\ = z\ + r, Z2 = Z2 + /3, -Z2 = Z2 + /3. 
Substituting these expressions in (4.5.1) one obtains 

fx(ii +r) = -zi, fi(z2 + '$) = -z 2 

from where 



-Zi - fir 



dz\ 
dt 

dz 2 ■■ 

fi— = -z 2 - (if) 



(36) 



(37) 



(38) 



So, for vector process (33) 



/ r (s,7,fi,Si,zi,z 2 ) = r, f l 3(s,-f,n,n,z 1 ,z 2 ) = 



(39) 



Expr 

[lOj . Let write the system (5) in the form 



Express values r and /3 (due to (5)) as continuous functions of state of the control object and regulator 

(40) 



Mr + Lf3cos/3 = Gu{t) + D(t) - Nr + L/3 2 sinf3 
Lcosj3f + 1(3 — Lgsin/3 — c(3 — nrcos/3. 

The solution of this system of linear algebraic equitions relative to variables r and (3 has a form 
I{Gu{t) - Br + Lsin(3(3 2 + D(t)) - Lcos(3(Lgsin(3(3 - Cf3 - brcosfJ) 



MI - L 2 cos(3 2 

M(Lgsinj3 - c/3 - bf cos/3) - Lcosf3(Gu(t) - Br + Lsinfi/3 + D(t)) 



(41) 



MI - L 2 cosf3 2 

In new variables the values r and (3 will be expressed with taking into account (10) as 

r = s — pO, (3 = H,cosf3 — (l 2 sinf3cos(3. 



(42) 
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In correspondence with (26) 



from where 



fl = difl - d 2 tg(3 - daQrsinfi + L 
(^ + as)(aL) 



yip — yip 
s = as 

ip z 



7 = -ay - AAu(t) - BD(t) 

-ay - AAu(t) - BD(t) + 2yL 2 sinf3 



(pL - 1)^ 



Here ip and its derivative are determined with using of relations 

ip = MI - L 2 cos(3 2 , ip = 2L 2 cosf3sin/3 = L 2 sin2(3 



(43) 



(44) 



(45) 



(46) 



and /3, in its turn is determined by relation (7). 

Now, turn out to verification of fulfilment of conditions of theorem 3 from To present the system 
considered in the form of 0, replace (38), (39) to system 



dzi 
~dt 



dz 2 
dt 



= ~Z2 - M2//3 



(47) 



So, system (35), (38) with functions 



fr = ? = s(s,y,Q, Q) - pQ(s,y,Cl,(l), 
fp = (3 = Q(s, 7, £1, (l)cosj3 — (l 2 sinj3cos(3 



(48) 



can be reduced to the form if denote x = col(s,y, z — col{z\, 



Z2 









X = 


1 
n 





( y/iP{(3) -as \ 
11 

-ay - J2 Xiuf 

i=l 

n = n 

\ n from (45) / 



g(x, z) = 



(49) 



z\ 



Tz 



-z\ 

— Z2 



h(x, z) 



s — pfl 

Qcosft — Qsin(3cos(3 



(50) 



At P2 — 0, z — it obtained the initial unperturbed system studied in As it was demonstrated in 
that work this system will have the dissipativity domain in the form of D — {(s, 7, O, f2), V(s , 7, fi, Cl) < 
C}, where V(s,y, tt, Cl) is the continuous piecewise positive definite in R m (in sense jl2|) function, 
satisfacting the condition a) of theorem 3 in 9]), i.e. 



dV(s,y,Cl,Cl) 



dt 



< — e, V(s, 7, fl, Cl) — > 00 at |s, 7, Cl, Cl\ — > 00, 



where S(R) is the bounded closed (hyper) surface: 

{(«,7,n,n), v(s,y,n,n) = R}, 



and R, e are the positive numbers. 
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Further, for system (38) 

dz 2 - (51) 

P~T~ = —z 2 — n(tlcos[3 — tl 2 sin(3cos(3) 

where s and are determined from (21), (25), it will take place a condition b) of the theorem 3 from 
0), and as W{z\, z 2 ) it can be taken quadratic form 

W=<z,Bz>, T T B + BY = —E (52) 

i.e. 

dW{z l Z2 ) <_jn N 2 ; ( 3 ^),(v Mie (o, Al ]) (W(z) ^oo \z\ ->oo). (53) 
at /ii 

Also from (32), (33) it will follow a fulfilment of condition c) of theorem 3 in Q. I.e. all conditions 
of theorem 3 will be hold and thus the system (5), (31)-(33), (38) will have the dissipativity domain 
D = {y,W p (x,z) < C}, independing of p G {Q,p] 3p > 0, W p (y) = V(s,7,fi,fi) + pW(zi,z 2 ) is the 
global Lyapunov function satisfacting due to the system considered the equality Q 



dW p (y) 



dt 



= -(\x\ 2 +p\z\ 2 ). 

M2=0 



Here y = (x,z), \y\ 2 = |(x, z)\ 2 — \x\ 2 + \z\ 2 , C p are some positive numbers. 
Then, it can be stated the following theorem. 

Theorem .Given system has attraction domain of equilibrium state x — 0, z — in the form W(x, z) < 
C at all < p < pi, and C — > oo pi — v 0. 

The proof is based on estimation of total time derivative of above introduced function 

w p (x,z) = w p (s,j,n,n;zi,z 2 ) = v p {s,-f,n,n) + pw(z 1 ,z 2 ) 

with using of inequality W p (x,z) < C as well as following from it inequality |(s,7, f2, f2)| < Rw p (C), 
where Rw (C) = sup\(s, 7, ft, Cl)\, W p (s, 7, fi, O; zi, z 2 ) < C following from conditions of the theorems 
estimations of right-hand sides 

J{s,^,Q.,tl;zi,z 2 ) < const(C)\(s,j,n,n)\\(zi,z 2 )\ 



h(8,j,Q,il;zi,z 2 ) < const(C)\(s,j,Cl,fl)\\(zi,z 2 )\ 

correct in this region. 



V. RESULTS OF NUMERICAL SIMULATION 

In work it was performed a numerical simulation for the case of closed-loop system with continuous 
control, described by equation (5). The behavior of mass-center cart coordinate as well as angular 
coordinate with increase of t evident about stabilization of the system for a short enough time interval 
at any values of parameters used. Moreover, with increase of parameter a oscillatory regime of decrease 
of angular as well as coordinate amplitude is changed to damping one. 

In present work, the numerical simulation of the behavior of the nonlinear system " inverted pendulum 
on a cart" with relay-type control under acting of continuously acting immeasurable perturbation was 
performed on the basis of MATLAB system. In particular, it was performed a numerical simulation 
of nonlinear dynamical system on the basis of the system (26) where as control force it was taken 
discontinuous relay-type (27) function with following parameter values a = 1; a = 0.5; k — 0.1; / = 
1; L = 1; M — 1.5; di = d 2 = 0*3 = 1; p m = 2. The initial values of variables were the following 
s = -0.7; 7 = 0.7; fl = 1.0; fl = 0.5. 
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For the case of absence of the control and external action (A = 0, B = 0) in the system all time depen- 
dencies for four variables s, 7, CI, Q are damping and asymptotically approaching zero but the character 
of approaching zero appears to be essentially different. The value s coming from the point with negative 
coordinate then becomes to be positive and approach zero from the positive value side. In contrast, 
the variable 7 is damping with increasing time up to zero being positive value beginning from t = 
(7(0) = 0.7). Then, time dependences fl(t) and Cl(t) are of oscillatory damping character. At t > 10 the 
system comes to equilibrium state s — 0, 7 = 0, O = 0, i7 = and rate of approaching the equilibrium 
state is determined by the friction force given by parameters L and k in the nonlinear system (35). 

The appearance of discontinuous relay-type control (27) in the system (A = 0.03) in absence of external 
action (B — 0) essentially decreases time of approach the equilibrium state in the system, moreover in 
7(i) dependence there appear horizontal regions 7 = corresponding to the sliding regime at the breaking 
surface in the state space. Note also that presence of relay-type control leads to significant decrease of 
time necessary for coming of the system in the equilibrium state fi = 0, r = 0. 

The numerical simulation in the case of external immeasurable sinusoidal perturbation in the system 
D(t) = sint in absence of the control (A = 0) all time dependences don't approach the stationary solution 
of the system (4.3.21) x = but are of oscillatory character with amplitude (Ax = 0.1) and period equal 
to that of external force D(t). The presence of discontinuous relay-type control (27) (A = 0.05) permit 
approach the stationary solution with sufficient accuracy (Ax — 0.03) already at t > 10. Such stabilization 
is determined, in essential degree, by appearance at j(t) dependence of horizontal regions corresponding 
to sliding regime at the breaking surface 7 = in the state space of the system (26) what don't permit 
to variables reach sufficient oscillation amplitude for time interval between horizontal regions. 

The presence of discontinuous relay-type control (27) in the case of action of external sinusoidal per- 
turbation leads to a fast stabilization of the system guiding the trajectory at the beaking surface 7 = of 
the system (26) and further it does not permit this trajectory essentially decline from this surface at pe- 
riodical sign change in external perturbation D(t). The essential fact here is that at more high-frequency 
perturbation the stabilization of the system occurs for a shorter time as compared with low-frequency 
(or even constant) external action. 



VI. CONCLUSION 



In this work it was performed a study of global asymptotic stability of nonlinear mechanical system 
"inverted pendulum on a controlled cart" in conditions of continuously acting perturbation with using 
method of two Lyapunov functions. It was demonstrated that at given choice of Lyapunov functions and 
comparison functions the conditions of the basic theorem of the method are fulfilled, i.e. system "inverted 
pendulum on a controlled cart" with discontinuous relay-type control, in conditions of continually acting 
perturbation is globally asymptotically stable. Here the solution is considered in the sense of Gelig et 
al. 8] . There is established the robustness property of above obtained stabilization algorithms relatively 
nonidcality of measurements of velocities of both cart and pendulum, since in practice there are no strict 
signal differentiators (instead of this the approximative velocity measurements are used). Also, using 
the method of two Lyapunov functions it is demonstrated that at given algorithm but with replace of 
ideal derivatives, the stability region persists, and its size at fi — > increases unlimitedly and in the limit 
coincides with the full phase space. 
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